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The concept of qudit (a d-level system) cluster state is proposed by generalizing the qubit cluster 
state (Phys. Rev. Lett. 86, 910 (2001)) according to the finite dimensional representations of 
quantum plane algebra. We demonstrate their quantum correlations and prove a theorem which 
guarantees the availability of the qudit cluster states in quantum computation. We explicitly con- 
struct the network to show the universality of the one-way computer based on the defined qudit 
cluster states and single-qudit measurement. And the corresponding protocol of implementing one- 
way quantum computer can be suggested with the high dimensional "Ising" model which can be 
found in many magnetic systems. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

Quantum computers can process computational tasks that are intractable with classical computers. The reason lies 
in the fact that quantum computing systems composed by qubits (two-level quantum systems) possess mysterious 
quantum coherence, such as entanglement (or quantum correlation], which has no counterpart in the classical realm 
Recently, an important kind of entangled states, cluster states 0, was introduced with remarkable property — the 
maximal connectedness, i.e., each pair of qubits can be projected onto maximally entanglement state with certainty 
by single-qubit measurements on all the other qubits. More surprisingly, it was shown that the cluster states can be 
used to build a one-way universal quantum computer, in which all the operations can be implemented by single-qubit 
measuments only As one of the key elements to realize such a scalable quantum computer, a celebrated theorem 
proved by R. Raussendorf and H. J. Briegel provides a simple criterion for the functioning of gate simulations on 
such quantum computer It was pointed out that the protocol of cluster state computers can be easily realized in 
practical physical systems since the creation of the cluster states needs only Ising- type interactions In fact, it has 
been demonstrated that the Hamiltonian with such interactions can be easily found in the solid state lattice systems 
with proper spin-spin interactions '5] , even in the system for cold atoms in optical lattices Q . 

Theoretically, it is natural to ask whether the concept of cluster states would have its counterpart in higher 
dimensional Hilbert space, i.e., the qudits cases, since most of the available physical systems can not be treated 
as a two-level systems even in an approximate way. The answer to this question is affirmative. Using a pair of 
non-commutative operators X and Z, which will be identified with the d dimensional irreducible representations of 
Manin's quantum plane algebra (QPA) in Sec. II Q, we find that the qudit cluster states \(f>)c can be defined as a 
common eigenstate of the tensor product operators 

Xl (g) Z,, (1) 

bGnbgh(a) 

where the lower indexes a and b denote qudit a and qudit b in the cluster, and index b is taken in the neighborhood 
of index a depending on the cluster structure. Based on this definition of qudit cluster states, we further show that 
we can build an one way universal quantum computer by explicitly demonstrating how to construct all single qudit 
unitary gates and one imprimitive two-qudit gate. 

We organize our paper as follows: We first briefly review the finite dimensional representations of quantum plane 
algebra in Sec. II, which provides the main mathematical tools in this article. As a non-trivial generalization of 
the qubit cluster state, the qudit cluster state is defined according to the quantum plane algebra, and its essential 
properties of quantum correlations are analyzed in Sec. III. In Sec. IV, the one-way quantum computer based on qudit 
cluster states is proposed with crucial supports from the proof of a central theorem. Similar to that for qubit case 
0, this theorem guarantees the functioning of gate simulations on those computers. Its proof depends on the subtle 
understanding about the quantum correlations of the qudit cluster states. In Sec. V, the universality of the qudit 
cluster state computers is proved through the explicit construction of all the single- and one two-qudit logic gates. 
Finally, we give our conclusion in the last section. 
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II. FINITE DIMENSIONAL REPRESENTAIONS OF QUANTUM PLANE ALGEBRA 



In this section, we will review the mathematics for finite dimensional representations of Manin's quantum plane alge- 
bra, which will be the main mathematical tool to describe not only qudit cluster states but also unitary transformation 
on the Hilbert space. The Manin's quantum plane is defined by 



XZ = qZX, 



(2) 



where g is a complex number. Mathematically, it can be proved that the associated algebra generated by Z, X 
possesses a d— dimensional irreducible representation only for q"^ — 1 In this article, we take q = qd = e'"r. This 
special case is first introduced by Weyl [13 and whose completeness is first proved by Schwinger . Obviously, when 
d = 1, q = 1, X and Z can be regarded as the ordinary coordinates of plane; When d = 2, q = —1, X and Z can 
be identified with the Pauli matrices ax and cr^. In this sense Z and X can be regarded as the so-called "generalized 
Pauh operators" 0, i, i, E lli IjJ ■ 

In fact, when q = qd, Z^ and X^ commutate with the algebra generators, so they belong to the center of QPA. The 
Shur's lemma tells us Z'^ and X'^ are constants multiples of the d— dimensional identity matrix, i.e., Z'^ = zl and X^ 
= xl. In general we can normalize them to the identity. Since the complex field C is algebraically closed, there must 
exist an eigen-state |0), which satisfies 



Z|0) = |0). 

According to Eq. 0, we obtain all the eigenvalue equations for operator Z 

Z\k)=ql\k), (keZd), 

where \k) — X'^''\0). This also impHes 

X\k) = \k-l). 



(3) 
(4) 
(5) 



In the Z-diagonal representation, the matrices of X and Z are: 



Z = 



1 

qd 












4-' 



(6) 



X 



1 
I 





1 



• 

• 



1 




(7) 



From Eq. lO, we have 



where 



x\xm = ii(o)), 



fc=0 



Similar to the eigenvalue equation of Z, we have 



where 



X\xm^q=Mj)) 



.(j))=Z^|x(0)) = ^^gf|fc). 



(8) 



(9) 



(10) 



(11) 
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Corresponding to this representation, we can also define Z{d) algebra as generated by Z, X. Its all basis elements 

B = {Z^X'',{j,keZd)}, (12) 
are called unitary operator bases in Ref. 'I4l. The general commutation relations for any two basis elements are 

X^Z'' ^ qfz'^X'. (13) 

In addition, we can replace the generators Z and X with two other elements in the basis. First, Let (m, n) be the 
greatest common factor of integers m and n. Then if (mi, ni) = 1, (rrii, ni e Zd), we can take 

^^^-^mi"i^mi^ni^ (14) 

where the factor before Z'^X'^ makes Z have the same eigenvalues with Z . To maintain Eq. (jSJ, we take 

X = (15) 

where (7712,^2) — 1, (w2,f72 G Zd)^ and min2 — m2ni — 1. From another viewpoint, Z and X defines a unitary 
transformation U 

Z = UZU\ X = C/XC/t. (16) 

By the above definition, it is easy to check that all this kind of unitary transformations form a group. In fact, it is 
the so-called Clifford group, which is a useful concept in universal quantum computation. 

III. QUDIT CLUSTER STATES IN QUANTUM PLANE 

To generalize the concept of qubit cluster states to qudit cases, we first restrict ourselves to one dimensional lattices 
for the sake of the conceptual simplicity. First, let us recall the definition of one-dimensional cluster states for N 
qubits. For a A^-site lattice, each qubit is attached to a site. As a novel multi-qubit entanglement state, the cluster 
state is written as 



1 A 



2N/2 yy 

a—1 



(|0)a-Ml)„(a,)„+i), (17) 



where {ai)a {i — x, y, z) are the Pauli matrices assigned for site a in the lattice, and 

a,|s) = (-l)^|s),(se{0,l}). 

Analogy to Eq. I|17(l . it is natural to conjuncture that the qudit cluster state in one dimension as 

N /d-i \ 

i<^)c = ^(8) Eifc)'^^'+i ' (18) 

0=1 \fc=0 / 

where 

Za\k)a=qi\k)a.'ia. (19) 

Now we present one of our main result: 

Theorem 1 The qudit cluster state in one dimension defined by Eq. ijjj))) is a common eigen-state, with eigen-values 
being equal to 1, of the operators X^ Zi, , i.e., 

6(znbgh(a) 

Xl (g) Z,\<j>)c^\4>)c, (20) 

6Gnbgh(a) 

where 

{2}, a=l, 
nbgh(a) ^{ {N~l}, a^N, (21) 

{a- l,a-f 1}, {l,N}. 
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Proof. To prove this theorem, we notice that qudit cluster state (|18|l can be constructed in the following procedure. 
Wc first prepare a product state 

N 

l+>=(R)l^(0))a, 



then apply a unitary transformation 



5= n ^-ft (22) 



fc-a=l 

to the state |+). Here Sab is defined by as an intertwining operator 

Sab\j)a\k)b^q^''\j)a\k)b. (23) 

It is easy to prove 

\^)c = S\+). (24) 

Since ^11+) = |+), it is easy to check 

SXlS^\<j,)c = \^)c- (25) 

The next step is to prove 

sxls^^xl (g) Zb 

To this end, for a, 6, c G {1, A^}, we observe that 

SabXlsl = ® Zb, (26) 

SabXlSl^Za^Xl (27) 

SabXlsl ^XlycilaM, (28) 

and 

SabZcSl, = Z,, Vc. (29) 

This completes the proof. □ 

Although the above proof is restricted to one dimensional cluster, it is convenient to generalize from one-dimensional 
qudit cluster to more complex clusters whether in two or three dimensional space. In fact, for a general cluster C 
with one qudit on each site, the cluster state \(j))c is defined by the following eigen-equations: 

Xl (g) Zb\(b)c^\^)c- (30) 

??Giibgh(a) 

Formally, the definition of a general cluster is the same as one dimensional case. Different clusters correspond to 
different relations of neighbours. 

Now we will discuss the properties of quantum correlations in the above cluster states under single qudit measure- 
ments. For simplicity, we still restrict ourselves to one dimensional case. 

First, let us discuss how to describe a von Neumann measurement for single qudit. Although the generator Z (or 
X) is not Hermitian, i.e., the eigenvalues of Z (or X) are not real, the non-degenerate eigenstates of Z (or X) can 
still form a complete orthogonal basis of Hd. Therefore they can be used to define a von Neumann measurement. For 
example, when we make a measurement marked by Z, we mean that we can obtain different results corresponding to 
different eigenstates of Z. 
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Next, we discuss the minimal number of single qudit measurements needed to destroy all the quantum correlations 
in qudit cluster states. For the one dimensional cluster state, we find all quantum entanglement will be destroyed by 
measuring Zja, (a = 1, 2, • • • , [7V/2]). 

Finally, we will discuss the most remarkable property - maximally connected - of qudit cluster state, i.e., each 
pair of qudits in the cluster can be projected onto maximally entanglement state with certainty by single-qudit 
measurement with all the other qudits. In fact, to project arbitrary two qudits in one dimensional cluster onto 
maximally entanglement state, we only need to measure X for the qudits between them, and Z for all the other 
qudits. It is easy to find that two or three dimensional qudit cluster states is also maximally connected. We only need 
to find a one-dimensional path connecting these two qudits and measure Z for all the other qudits which are not on 
the path, which reduces the two or three dimensional problem to one dimension. 

We now discuss the problem of physical implementation of our one-way quantum computer. Physically the qudit 
cluster state (fl^ can be created by Hamiltonian 



H=-hgY,Ni^^Nlf\{9>^) (31) 

(a,6) 

where (a, b) denotes sites a and b are nearest neighbors in the cluster; and N'^^'^ is defined as 

iV(^) =^|fc)fc(fc| (32) 

fe=0 

Then we find the intertwining operator S in Eq. H22|l has the the explicit form 

S = exp ^— 

where the evolution time tc = 



S = expl--Htc] , (33) 



To associate with more familiar Hamitonian in physics, let us define the spin--^^^— !^ operator of z direction 



= iv(^) - (34) 



Then we can rewrite Eq. H31|l as 

a (a.b) 

where Va is the number of nearest neighbors for qudit a in the cluster. Obviously, the interaction Hamiltonian 

Hi = -hg'^{sz)a{sz)b, (36) 

{a,b) 

which is the ferromagnetic Ising type interaction with spin-^^. 

Those kinds of Ising model, other than the usual spin-1/2 Ising model, has been one of the most actively studied 
systems in condenced matter and statistical physics due to their rich variety of critical and multicritical phenomena. 
For example, the spin - 1 Ising model with nearest- neighbor interactions and a single-ion potential is know as the 
Blume-Emery-Griffiths (BEG) model [13, the spin - 3/2 Ising model was introduced to explain phase transitions in 
DyVOA and its phase diagrams were obtained within the mean-field approximation |l8j . In addition, higher spin 
Ising models can be associated with the magnetic properties of artificially fabricated superlattices. Such lattices 
consist of two or more ferromagnetic materials have been widely studied over the years, because their physical 
properties differ dramatically from simple solids formed from the same materials. The development of film deposition 
techniques has aroused great interest in the synthesis and study of superlattices in other materials. A number of 
exDerimental tia..20.. .21. 22.J and theoretical works |23^ .24. ■25„ ,.26„ .27, ■2& . .29, .30j have been devoted to those directions. 
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IV. MEASUREMENT BASED QUANTUM COMPUTATION WITH QUDIT 

As discussed above, the qudit cluster states exhibit the same features in quantum entanglement as that for the 
qubit cluster states. Then a question arises naturally: Can these natures of quantum correlations be available for 
constructing the universal quantum computations? We will give an affirmative answer to this question in the following 
two sections. 

In this section, we further generalize the basic concept of "single qudit quantum measurement" and the correspond- 
ing measurement based quantum computation (MBQC) on qubit clusters to that on qudit clusters. Along the line to 
construct MBQC for the qubit case, we will formulate the corresponding theorem which relates unitary transformation 
to quantum entanglement exhibited by the qudit cluster states. Quantum computations with qudit clusters inherit 
all basic concepts of those with qubit clusters. They include the basic procedure of simulation of any unitary gate, 
the concatenation of gate simulation and the method to deal with the random measurement results. Here we will give 
a c?-dimensional parallel theorem as generalization of the central theorem 1 in Ref. Q . 

Before formulating our central theorem, let us introduce the basic elements for quantum computing with qudit 
clusters. The main task of quantum computing with qudits is to simulate arbitrary quantum gate g defined on n 
qudits Hilbert space. For this purpose, the first step for quantum computing on qudit clusters is to find out a proper 
cluster C{g). Then we divide it into three sub-clusters: the input cluster Ci{g), the body cluster CM{g), and the output 
cluster Co{g)- As usual we require that the input and output clusters have the same rank (i.e. the same number of 
qudits), |C/((7)| — \Co{g)\ = n. Then we prepare the initial state in 

(37) 

In this first step of quantum computing, we entangle the qudits on the qudit cluster by using the cluster state 
generator S, i.e., 

|$(in))c(,) -5|f(in))c(,). (38) 

This step brings the structure information of the qudit cluster into our computing process, and thus relates it with 
the corresponding qudit cluster state. 

The second step is to measure all qudits on the cluster in special space-time dependent basis according to a given 
measurement pattern (MP). The definition of MP is given as follows. 

Definition 1 A measurement pattern M.c on a cluster C is a set of unitary matrixes 

Mc = {uaZaui\aeC,UaeSU{d)}, (39) 

which determines the one-qudit measured operators A^i"' on C, with the explicit form 

5=0 

If this measurement pattern Mc operates on the initial state |$(in))c(g), the set of measurement outcomes 

{■s}c^{saeZd\aeC} (41) 
is obtained. Then, modulo norm factor, the resulting state |5'a4)c is given by 

\'^%\)-p1Z\'^(^^))c(,), (42) 

where the pure state projection 

PMc=^^k\s)kk{s\ul (43) 
feec 

It is worthy pointing out that we always measure X for the input qudits and Z for the output, which is independent 
of gate g, i.e., 

Mcig) = {X,,ieCi{g)}, (44) 

Mcoig) = {Z^,ieCo{g)}. (45) 

Then we reach the final step, to associate the measurement values with the result of gate g acting on the initial state. 

From the above standard procedure of qudit clusters quantum computation, we learn that what is crucial for this 
scheme is to associate a given gate g with a measurement pattern. Although by now we have no general optimal 
operational procedure to do this for practical problem, the following theorem provides a useful tool in realizing specific 
gates on the qudit clusters. 
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Theorem 2 Suppose that the state |V')c(g) = (g) \^)c{g) obeys the 2n eigenvalue equations 

^c.te)4f^^^f^0cote) l^)c(g) = (46) 
Z^c,ig)AUZ^U^),^^^^Wcig) = 9,'^"IV')c(,), (47) 
with Xx,i, Xz.i G Zd and 1 < i < n. Then, according to the above standard quantum computing procedure, we have 

PmI^^^PmI,,,, l*(i^)>c(.) « (ni^»))cH.)uc„te)lV'(ont))c,(,), (48) 
where the input and output state in the simulation of g are related via 

|i^(out)) = C/C/s|V'(in)), (49) 

where Uj] is a byproduct operator given by 

n 

Us^ (g) (Z,)-^-'-^'(X,)^-. (50) 
Proof. Let us begin with the case when 

min))c,(g) = \{t})cAg), (51) 

with 

{t} = {ilia •••*«}. 
To associate with |V')c(s)i the initial input state is written as 

|V'(in))cH,) = (Vd)"Fj^! |+)c,(.), (52) 

Cj (s) 

where 

M'c,^g^^{Z.,ieCi{g)}. (53) 

In terms of Eqs. we have 

i{g)UCM{g) 

To find out the equations for the final state \->p{ou.t)) i^^^^^y we take PjCi], ^ ^-^m' acting on both sides of Eq. ijl^ 
and Eq. (|T7|l : 

(t^^*t^^)co(.) l^(°^^))co(.) - ir''^''\^i°^^))coigV (55) 

(U^^U^)coig) l^(°^^))co(.) - lT'"\^i°^^))coigV (56) 
where the input state for |?/'(out)) 

min))=Xj\{t}). (57) 

Before drawing a conclusion, we need to check the final state is not a zero vector. In fact, from Eqs. H4t)|) and 
(|47ll . the state U''\ip)c{g) is the simutaneous eigen-state of operators ^C/(s),i-^Co(s),j ^^'^ Z'^ Ci{g),iZco(g),i- We can 
directly evaluate this state and find out that it has every components in Z-diagonal representation of the input part. 
Consequently the final state is indeed a nonzero vector. So from Eq. H55|l and Eq. H56|l. we obtain 

|V(out))c„(3) = e^^Wt/t/s \{t})co(g)- (58) 
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To further determine the relation between the output state and the input state, let us consider the other case when 

|V'(in))cH,) = l+)c.(9)- (59) 

At this time, the final state 

\{^i})ci(g)uCM{g) I ^' (out)) Co (g) = ^it)(^)l^)c(.)- (60) 
Let f J ^ apply on both sides of Eq. H46() and Eq. (|47|l , we obtain 

(^^^U^)co(,) \^'io^^))coi9) = qr''"Wiont))coi,^, (61) 

i^^^^^)coi,) \^'i°^^))co(9) = q^'-'\^y (out)) co(,), (62) 
where the input state for |?/''(out)) 

(in)) =zi\+). (63) 

From the above equations, we obtain 

|V'(out))co(,) =e^'^C/t/s|+)c„(g). (64) 

Substitute Eq. ^ into Eq. 

|V'(out))co(,) = C/C/Ed-^ E \{t})coi9)- (65) 

{t} 

Comparing Eq. H64|l to Eq. (|55|) . we obtain 

e^''^*) = e*x. (66) 

This completes the proof. □ 

This theorem tells us that, since the cluster states have remarkable quantum correlations, which has been discussed 
in Sec. II, they play an essential role in the realization of arbitrary unitary gates. More precisely, as long as one 
cluster can be used to process a unitary gate for the cluster state, it will work for arbitrary input states. Therefore, it 
is sufficient to check the conditions for the cluster states, i.e., Eq. H46(l and Eq. H47|l. To be emphasized, the special 
features of the above theorem arising from qudits is that it is expressed not only in terms of unitary operators X and 
Z , but also in terms of their conjugates. For d = 2, it exactly reduces to the Theorem 1 of Ref. [J. 

Before using the theorem to construct a specific unitary gate, we need to explain how to deal with the byproduct 
part Us- The basic idea is to move J/s to the front of U according to the commutation relations between and U. 
To complete this operation, the general strategy is to divide measurements into several steps, in which the subsequent 
measurements depend on the results of previous measurements. In the next section, we will use specific examples to 
demonstrate how to construct all basic elementary gates with the help of this theorem. 



V. UNIVERSALITY OF QUDIT CLUSTER QUANTUM COMPUTATION 

It is well known that, for qubit quantum computing network, a finite collection of one qubit unitary operations and 
CNOT gate is enough to construct any unitary transformation in the network. This conclusion has been shown to 
work well for qudit case |l6j | . More precisely, the collection of all one-qudit gates and any one imprimitive two-qudit 
gate is exactly universal for arbitrary quantum computing, where a primitive two-qudit gate is such a gate which 
maps all separate states to separate states. Therefore, in order to prove universality of quantum computation with 
qudit clusters, we only need to construct those basic element gates, and then to integrate them to realize arbitrary 
unitary gate. In this section, based on Theorem[5]in Sec. Ill, we will build certain qudit clusters to simulate the basic 
elementary gates, i.e., any one-qudit unitary operation and one imprimitive two-qudit unitary transformation. 
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A. Realizations of single qudit unitary transformations 

Let us start with single qudit gates. First, we introduce a proposition for any unitary transformation in d dimensional 
Hilbert space. 

Proposition 1 Let {Ni, (i G Z^2_i)} be a Hermitian basis of the operator space for d dimensional Hilbert space, 
then any unitary transformation U has the form 

U = q,^^'^^^^=l[q/^^% (67) 
i 

where Ui and f3i are real numbers. 

The first equality is obvious, and we refer to Ref. for the proof of the second equality. 

By the above proposition, we can divide any qudit gate into a product of more simple basic ones, i.e., single 
parameter unitary transformations. Now, we need to find d x d independent Ni to simulate all unitary gates for one 
qudit. To use the definition of qudit cluster states, we expect that the single parameter unitary transformations must 
have deep relations with the basis elements of QPA. In fact, we find such a way to introduce d x d one parameter 
unitary transformations. This originates from that the fact that some of the basis elements of QPA can be used to 
define a state basis. We find that all the unitary transformations that do not change every basis state up to a phase 
are defined by the property of multi- values for complex functions. For example, for operator Z, we can define 

^''({w)='?^™^(/3e7^), (68) 

where 

d-l 

N{Z,{m}) = |n)(n + m„d)(n|, (Vm„ £ Z). (69) 

Although the above definitions include infinity unitary transformations, there are only d independent ones, which can 
be used to describe the following type of unitary transformations 

Uz{{a})\n) = (Vn G Zd,a G 7^). (70) 

Obviously, these d independent unitary transformations, also the corresponding N{Z,{m}), can take the place of 
{Z'^,n G d} in the unitary basis. 

A similar argument can be generalized to Z, which is defined in Sec. II. More precisely, for operator Z, we define 

^^(W) = 9r'^^^"", (/3e7e), (71) 

where 

d-l 

Af(Z,{m}) = ^|n(Z))(n + m„d)(n(Z)|, (Vm„ G Z), (72) 

n=0 

with 

Z\n{Z))^q2\n{Z)). (73) 

In the following, we will show that we can select d independent Hermitian operators from N{Z, {m}). 
When d is a prime number, a convenient choice is to take Z from the operator set 

{Z,X,ZX,ZX^,---,ZX'^-^}. (74) 

Because each Z defines d—1 independent N{Z, {m}) besides the identity, we obtain d^ (= {d—l){d+l) + l) independent 
Hermitian operators. 

When d is not a prime number, we can choose the independent N{Z, {m}) by the following procedure. First, we 
take Z from {Z, X, ZX}, and we can obtain 3{d — 1) independent N{Z , {m}) besides the identity, which can take 
the place of the set of basis elements 



S = {Z", X", {n G Zd)). 



(75) 
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Then we take an element Z ^ S, find out the elements in {Z", n € Zd\ that is not in the set S, add these elements 
into S, and take the new independent N{Z , {to}), whose number is the number of new elements in set S. We repeat 
the above step until set S — B, then we obtain cP independent Herniitian operators. 

Therefore if we can do all the above basic unitary transformations Z^{{m}, we will declare that we can do all single 
qudit unitary transformations. Here we adopt the following strategy: First, we realize Z"({to} on a five qudit cluster, 
as a basic single qudit transformation. To associate with all the other single qudit unitary transforamtion Z'^^y, we 
observe an useful fact 

Z"({to}) = C/Z"({TO})^7^ (76) 

where U (or W) can be taken as an element in the Clifford group as defined in Sec. II. Based on this observation, 
we further discuss how to implement these elements in the Clifford group on clusters. Then we only need to connect 
these clusters in the given order to realized unitary gate Z°'({m}). If we succeed to pass the above procedure, in 
principle we can make any single qudit unitary gate in SU{d). 

Now, let us realize these basic unitary transformations with a specifically-designed qudit cluster, which will be given 
as follows. 



1. Five qudit cluster realization of X" {{m}) and Z"{{m}) 



In this subsection, we aim to realize the basic single qudit unitary transformation Z"({to}) on a five qudit cluster 
designed as in Fig. ^ which is a linear array of five qudits. By the way, we also use the same cluster to implement 
X"({to}), which demonstrates that the same cluster with different measurement patterns can realize different unitary 
transformations. The corresponding cluster state is defined by the following system of equations: 

h Oit 



FIG. 1: Five qudit cluster used in realization of X°'{{m}) and Z°'{{m}). A circle represents one qudit, number n in the circle 
means the n-th qudit, in or out denote the input or output part of the cluster, and two qudits which are connected by a line 
is neighbors. 



It follows from Eqs. (|77I81|) that 



Form Eq. (|HTT) . we obtain 





^)c = 1 


0)c 


ZiXlz^lc, 


^)c = 1 


0)c 


Z2XlZi\c, 


^)c = 1 


4>)c 


Z^X\Zt^\c, 


^)c = 1 


<P)c 


ZiXl\c, 


^)c = 1 


4>)c 



XiX'lX^] 

z\x2x\z<^\ 



zt{{m})X^{{m}mc 



(77) 
(78) 
(79) 
(80) 
(81) 



(82) 
(83) 



(84) 



Here we emphasizes that when we obtain the above equation, we have used the following condition On {to}. If 
n4 + rig = {Mod{d)), then 



rii + TO„4 +^5-1- TO„5 = 0, (n4, n5 G Zd). 



(85) 
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From the above four equations, we have 

X,XI (X^ {{m})X,xl" ({m})) \<t>)c = \<t>)c. (86) 
Zlx, [zl-{{m})xlZ2{{m})) (x^{{m})Z,xl-{{m})) |0)c = |</.)c. (87) 

For a measurement pattern {X2, Xj, z]^{{m\)x\Z2{{m\)}^ Theorem |5] concludes that the simulated unitary trans- 
formation is X^{{m\)UY., where 

Because C/s depends on the measurement results and can not be moved to the front of X^{{m}) trivially, different 
measurement results lead to different unitary transformations. In order to realize unique gate X^({m}), we complete 
the measurement in two steps: In the first step, we measure {^Xi^X2tX^}. Then, when the obtained values are si , 
(= S3) and S2, the byproduct operator reads C/s — Z~^^~''^ X^' . At this time, \z is still unknown since it depends 
on the measured result S4. However, notice that Z'!^^^^'^'^ X^{{'m\)Zt,^'^^^ is diagonal in the X representation, we have 

Z-'^-^^X5"({m})Zf = n ^"''"({W)- (88) 

{m} 

From Eq. (|88() . we can obtain d equations of {a^J^j}, which determine the values of 

Now we make a new choice depending on the known measurement results. Also from Eqs. (|77I81|) . we obtain the 
following equation instead of Eq. 18 7|) 

Z\X, ( n zf'-h{^})xl n ({m}) 

\{m} {m} 

n {{m})Z, n xf'"-' ({m}) I |<^)c - |0)c. (89) 

Jm} {m} 

Measuring the fourth qudit in the base 



SI S3 



{m} {m} 



we obtain the value S4 and therefore \z =52+54. According to Theorem |5J we obtain the final operation 

'■5 ^5 



UU^ = Z-'''-''Xl'+''X^{{m}). (90) 



Finally, by measuring Z5 we obtain the correct result 

S = S5 + S2 + S4- (91) 

The above equation does not mean that the final result depends only on the measurements to the second, fourth, 
and fifth qudit, because different values of si and S3 correspond to different measurement on the fourth qudit. 

Based on the same cluster, we can also implement the single qudit rotation Z"({m}). Similarly, we first measure 
{Xi,X2,xl}- Then from Eqs. (|77I81|) . we obtain 



n n ^5 ({"^}) = \<i>)c, (92) 

\ m m / 

zlX2Xlz^\cf,)c = |0)c. (93) 

where {q;|^|} is determined by 

X^,^+^-Z^{{m})Xr^-^- = n zf-^ ({m}). (94) 

{m} 
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At the same time we make another measurement on 

= 2=4 



According to theorem [3 we conclude that the simulated unitary transformation indeed is 
The correct result and the measurement values are also related by Eq. H91|l . 



2. Realizations of single qudit elements in Clifford group 

As implied in Eq. H76|l . we only need to realize the single qudit elements in the Clifford group. It is easy to show 
that not all elements in the Clifford group are required. In fact, we only need do the elements defined as 



where 



with 



Id 



= z. 



(m, n) = 1. 



(95) 
(96) 
(97) 



We will show that we can do all the above Clifford unitary transformations through a series of four basic types of 
unitary transformations. The first is defined as 



The second is defined as 



The third is defined as 



The last is defined as 



VZV'' 
VXV'' 

WZW'^ 

wxw^ 



Id 

X; 



ZX"". 



<ld ' zx, 

X: 



= z, 

= Qd ' ZX. 



(99) 

(100) 
(101) 

(102) 
(103) 

(104) 
(105) 



Theorem 3 For any unitary transformation JJ"^" partially defined by Eq. it can be factorized into a product of 

a series of the above four basic unitary transformations. 

Proof. Let us prove it by induction. We denote Nd = {1, 2, • • • , d}. When m = 1 or n = 1, it is the first or the 
second types of unitary transformations. Let us suppose the above theorem is valid at m G Nd or n G Nd, i.e., we 
can do 



f/™" {moine Nd, im,n) = 1). 
For arbitrary positive integer n that is relatively prime to d + 1, there exists 

n^i{d+l) + n', {i e Z^, n' e Nd). 



(106) 
(107) 
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Obviously, 

n' € Nd, {d+l,n') = 1. (108) 
According to the assumption, we can do f/(''+i)" . Then let V operate i times, 

YiUid-+l)n'Y^!^ ^ jj{d+l)n^ (109) 

According to the assumption, we can also do [/" Then let W operate i times, 

^ij|n'(<i+l)|yt' ^ (110) 

Therefore the theorem is valid for m or n E N^+i- This completes the proof. □ 

Now we come to constructing the four basic unitary transformations. We will prove that the first (including the 
third) and the fourth can be realized on the five qudit cluster as Fig. ^ According to Eqs. H77I81|I . we have 

X,xlX5\<f>)c = |0)c, (111) 

ZtX2(Z4"X4)tZ5X5"|0)c = |0)c. (112) 

According to theorem|2 when the measurement pattern is {X2, xl, [qd '^'^ Z^^ Xn)'^}, the corresponding unitary 
transformation is 

q^^^{si+S3)n ^-si-S3 j^S2+Si-n(si+S3)jjln (113) 

Also from Eqs. (|77I81|I . we obtain 

XiZ:,xlxlz^X^\<t>)c = \4>)c, (114) 
zlX2Xlz^\ct,)c = |0)c. (115) 

According to theorem|21 when the measurement pattern is {X2^ qd^~^Z:iXl, xl}, the corresponding unitary trans- 
formation is 

To realize the unitary transformation [/"^, we need a cluster composed of six qudits as Fig. |21 The cluster state is 



fi Oit 



FIG. 2: Six qudit cluster used in realization of U"^ . Meanings of the symbols in this Figure are the same as Fig. 
defined by the following system of equations 



It follows from the above equations that 



xlZ2\^c) - 



ZiX^Zsl 

Z2XlZ4\ 

ZzX\z<^\ 

z^xlz^\ 



XiXlx<:,zl\ 
z\X2Z4'^ x\xl Zq^XqI 



When the measurement pattern is {X2, X^, qd 2 "^4"^!, X5}, the corresponding unitary transformation is 



(117) 
(118) 
(119) 
(120) 
(121) 
(122) 

(123) 
(124) 



(si+S3)(s2+S4-nS5 + - 



7n(s5 — Si— S3) — S2 — S4 j^ — Si — S3 j-fnl 



(125) 
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B. Realization of an imprimitive two qudit gate 



Now we come to the construction for simulating two qudit operations. The cluster composed of six qudits as Fig. 








Oct 

6 




FIG. 3: Six qudit cluster used in realization of an imprimitive two qudit gate T. Meanings of the symbols in this Figure are 
the same as Fig. 

01 is considered with the following system of equations 





be) = 1 




(126) 


^2 ^4 1 


^>c = 1 


0)c, 


(127) 


Zi z^ 1 


^>c = 1 


0)c, 


(128) 


Z2ZzX\Zq\ 


^>c = 1 


4')c, 


(129) 


ZaXll 


^>c = 1 




(130) 


Z4,X\\ 


^>c = 1 


0>C. 


(131) 



It follows from the above equations that 



X2Xl\ 



^i^aZlx, 

zlx^zlxl\4>)c 



(132) 
(133) 
(134) 
(135) 



By measuring the system according to the measurement pattern the simulated two qudit gate 

T satisfies, and is also defined by 



TXsTt = Xl 
TXeT^ = xl 
TZ,T^ = Zlxl 

tZqT^ = zlxl 

According to Theorem 12 the above measurement pattern realize the following unitary gate 



rpy — Sl VS3 y — S2 VS4 



(136) 
(137) 
(138) 
(139) 

(140) 



The next task is to prove that T is an imprimitive two qudit operation. Ref. 16] tells us that a two gate V 
is primitive if and only if V = Si ^ S2 or V — (Si S2)P- Here, and ^2 are different single qudit operators, 
P is the interchanging operator obeying P\x) (8) \y) — \y) <Xi \x). Based on this fact, we can easily conclude that a 
primitive operator always maps a single qudit operator to another single qudit operator. Obviously, the above two 
qudit operator T is imprimitive. Another way to prove it is to evaluate the unitary transformation directly. Then we 
can find that it maps all Z (g) Z— bases to the maximum entangle states. 

As demonstrated in this section, any single qudit unitary gate and one imprimitive two qudit gate can be realized 
on qudit clusters. Therefore, the measurement based quantum computing on qudit clusters is universal. 
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VI. CONCLUSIONS 

We have introduced the concept of qudit cluster state in terms of finite dimensional representations of QPA. Based 
on these qudit cluster states, we have built all the elements of qudit clusters needed for implementation of universal 
measurement-based quantum computations. With generalizations of cluster states and measurement patterns, most 
of the results in qubit cluster can work well for qudit clusters in parallel ways. We also show that there still exists 
the celebrated theorem guaranteeing the availability of qudit cluster states to quantum computations. To prove the 
universality of this quantum computation, we show that we can implement all single qudit unitary transformations 
and one imprimitive two qudit gate on specific qudit clusters. In addition, we propose to build a one-way universal 
quantum computer with qudit cluster states practically since the high dimensional "Ising" model can be used to 
generate such cluster state dynamically. 
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